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ABSTRACT 

We use about 200,000 FGK type main-sequence stars from the LAMOST DR1 data to map the 
local stellar kinematics. With the velocity de-projection technique, we are able to derive the averaged 
3 dimensional velocity and velocity ellipsoids using only the line-of-sight velocity for the stars with 
various effective temperatures within 100 < \z\ < 500 pc. Using the mean velocities of the cool 
stars, we derive the solar motion of (U®, V®, Hq)=( 9.58±2.39, 10.52±1.96, 7.01±1.67) kms -1 with 
respect to the local standard of rest. Moreover, we find that the stars with T e g> 6000 K show a 
net asymmetric motion of ~ 3 kms -1 in (W) compared to the stars with T e g< 6000K. And their 
azimuthal velocity increases when \z\ increases. This peculiar motion in the warmer stars is likely 
because they are young and not completely relaxed, although other reasons, such as the resonance 
induced by the central rotating bar or the spiral structures, and the perturbation of 
the merging dwarf galaxies, can not be ruled out. The derived velocity dispersions and cross 
terms for the data are approximately consistent with previous studies. We also find that the vertical 
gradients of ajj and ay are larger than that of aw . And the vertical gradient of au shows clear 
correlation with T e g, while the other two do not. Finally, our sample shows vertex deviation of about 
11°, at 300 < \z\ < 500pc, but roughly zero at 100 < \z\ < 300pc. 

Subject headings: Galaxy: disk Galaxy: evolution - Galaxy: fundamental parameters - Galaxy: 
kinematics and dynamics - solar neighborhood 


1. INTRODUCTION 

The velocity distribution of the stars in the solar neigh¬ 
borhood plays a key role in understanding the global 
structure, dynamical features, and the evolution of the 
Milky Way. Although it is often approximated with a 
multi-dimensional Gaussian profile, the velocity distri¬ 
bution of the stars in the solar neighborhood is actually 
very complicated. The mean value of the velocity dis¬ 
tribution should be around zero, given that the Galactic 
disk is in static state. Howev er, observat ions have found 
many substructures (iDehnen H998tlZhao. Zhao fo Chen 


2009: Sicbert et al. ll201lHAntoia et al. 2012: XiaetaL 


201411 , which may be associated with the perturbation of 
the Galactic bar and spiral arms, or belong to old tidal 
debri s of di s rupte d clusters or dwarf ga l axies (IDehnen 


2000; IFux I 1 20011 iFamaev et al. I 120051 lAntoia et al. 


20111 ). in the velocity distribution. These substructures 
may shift the mean velocity slightly away from zero by a 

few kms -1 ._ 

Recently, IWidrow et al. I (120121 1 found that the stel¬ 
lar number density is not symmetric with respect to the 
Galactic mid-plane, imply ing a vertical wave in the stel¬ 
lar disk. iGomez et al. 1 (1201311 inferred that the merg¬ 
ing event of the Sagittarius dwarf galaxy can induce such 
density waves. The vertical asymmetry in stellar count 
may be associated with vert ical asymmetry in velocity. 
Indeed, I Carlin et al. I ()2013H found from the LAMOST 
DR1 data that not only radial velocity but also the verti¬ 
cal velocity of the nearby stars are not symmetric. These 


new observational evidences challenge the current dy¬ 
namical models of the disk. At least some of the nearby 
disk stars are not in equilibrium, although the majority 
must be in the static state so that the Galactic disk can 
survive over billions of years. More investigation is re¬ 
quired to figure out how many and which types of stars 
contribute to the asymmetric motion, which is our main 
motivation of this work. 

The velocity distribution can be characterized by 
the velocity ellipsoid, which reflects the mass dis¬ 
tribution and evolution of the Milky Way, assum¬ 
ing that most of the detected stars are in equilib¬ 
rium. The earlie st st udy of_the atellar velocity ellip¬ 
soid was done bv iSchwarzschildl (119081) . From then on, 
many works have found that the age of stars is corre¬ 
lated with the velocity distribution. Specifically, older 
stars show larger velocity dispersion, and vice versa 
(iParenae ro Ill95flt lRoma.nl I* 9511119521 IDehnen fc Bin nevl 

19981 iQuillen fc Garnett l~ 2001b iNordstrom et al.l 12004 

Holmberg, Nordstrom fe Anderson 2Q0T etc.). This is 

usually thought to be because scattering of the disk 
stars increases over time. For the stars younger than 
~8 Gyr, the scattering is most likely due to encoun¬ 
ters with substr uctu res, e.g., g iant m olecu lar clouds, in 
the disk (jHolmberg. Nordstrom fc Andersonl2007l ). The 
age-velocity dispersion relation (AVR) reflects the evo¬ 
lution history of the Galactic disk. 

In an ideal axisymmetric disk, the velocity ellipsoid 
near the disk mid-plane should be aligned with the cylin- 
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drical coordinates. However, since the Galactic disk con¬ 
tains a rotational central bar and a number of spiral 
arms, the velocity ellipsoid will not match this ideal case. 
On the contrary, it deviates from the cylindrical coordi¬ 
nates near the Galactic mid-plane. iDehnen fc Binnevl 
( 1998 1 hereafter DB98) reported that the vertex devi¬ 
ation of the velocity ellipsoid measured from Hipparcos 
proper motion data is ~ 10°. iSrnith. Whiteoak fe Evaiisl 
(|2012L hereafter S12) confirmed that the vertex devia¬ 
tion of the metal-rich stars in the SDSS sample is con¬ 
sistent with DB98. Moreover, the authors also showed 
that the tilt angle relative to the Galactic mid-plane 
is about —10 to —15° for the velocity ellipsoid at z > 
0.5 kpc. Recently, lots of new works have been done with 
the data from larg e ste llar spectroscop i c sur veys, e.g., 
IBinnev et al.l (12014 ) and Sharma et al. I (120 Ml) ba sed on 


RAVE ( Kordopatis et al. 201 31) . IBovv et all (120 Ml) based 
on APOGEE ( Ahn et al. 112014) . and Recio-Blanco et ahl 
(|2014D based on Gaia-ESO survey dGilmore et al.ll20 1 2li . 

Since 2011, the LAMOST surv ey dCui et al.l l2012t 
iZhao et ahll2012t iDeng et al. II2012T) has been operating, 
and about 1 million stellar spectra with stellar param¬ 
eters have been released in the DR1 catalog. A large 
fraction of these stars are located within 1 kpc around 
the Sun, providing a vast resource to reveal details of 
the local stellar velocity distribution and give constraints 
on the dynamical structure and evolution of the Galactic 
disk. 

In this paper, we use the FGK type main-sequence 
stars selected from the LAMOST DR1 catalog to study 
the local velocity distribution. We use the effective tem¬ 
perature of these stars as a proxy for age and investigate 
the variation of the velocity distribution, including the 
mean velocity and the velocity ellipsoid, as a function 
of the effective temperature and height above/below the 
Galactic mid-plane. 

The paper is organized as follows. In Section 2, we 
describe how we select the samples. In Section 3, we in¬ 
troduce a velocity de-projection method to reconstruct 
the velocity and velocity ellipsoid in three dimensions 
from only the line-of-sight velocity. The method is then 
validated with GCS data and a mock catalog. The sys¬ 
tematic biases due to the uneven spatial sampling and 
the spatial variation of the velocity ellipsoid are then 
calibrated. The mean velocities and velocity ellipsoids in 
various effective temperature bins at different height is 
then derived and discussed in Sections 4 and 5. Finally, 
brief conclusions are drawn in Section 6. 


2. THE DATA 

2.1. The LAMOST Survey 

LAMOST is a quasi-meridian reflecting Schmidt tele¬ 
scope with an effective aperture of about 4m and a field 
of view of 5°, operated by the National Astronomical Ob¬ 
servatories, Chinese Academy of Sciences. It is a power¬ 
ful instrument to survey the sky with th e capab ility of 
recording 4000 s pectra simultaneously (ICui et al.l 120121 : 
IZhao et ahll2012f) . In its five-year survey plan, it will ob¬ 
tain a few millions stellar spectra in ab out 2 0 ,000 squar e 
degrees in the northern hemisphere (IDeng et al. Il2012il . 

The LAMOST Survey has internally delivered the first 
data release (DR1), which contains 2,204,860 spectra 
with a resolution of ~1800 covering a wavelength range 
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Fig. 1. Fractional distance residuals resulting from run¬ 
ning the ICarlin et al. I ((2015 ) distance code on a realiza¬ 
tions of the Besancon model in a field at ( l,b ) = (180°, 60°). 

of 3800 < A < 9100A. The catalog contains 1,085,404 
stellar spectra with estimated stellar atmospheric param¬ 
eters as well as line-of-sight velocities. The distances 
to stars with estimated effective temperature, surface 
gravity, and metalli city a re determined from isochrone 
fitting by ICarlin et al. I (I2015H . They developed two 
techniques, chi-square and Bayesian to estimate the ab¬ 
solute magnitude. The chi-square technique finds the 
best fit absolute magnitude which reaches the minimum 
of the chi-square between the measured T e g logy, 
[Fe/H] and the isochrone ones. The Bayesian tech¬ 
nique, on the other hand, considers both the selection 
effect and the different stellar populations along differ¬ 
ent line of sight. They verified that the accuracies of 
the distance estimates derived from both techniques are 
~ 20%. They also found that their method may under¬ 
estimate the distance of the giant stars by about 20%, 
but for main-sequence stars with 4000 <T e g< 7000 K, 
no correlation is found with T e g, as presented in 
Figure [l] Compared with the Hipparcos data, they ob¬ 
tained fairly good distance estimates with scatter of only 
~ 17%. Considering that only the FGK main-sequence 
stars are used in this work, the accuracy of the distance 
is sufficient for the studies of the local stellar kinematics 
within about 1 kpc. For such nearby stars, not significant 
difference is found between the two techniques. There¬ 
fore, we adopt the distance estimates from the chi-square 
technique in this work. 

2.2. The sample selection 

We select nearby FGK type main-sequence stars to in¬ 
vestigate the kinematics of the solar neighborhood. In or¬ 
der to select main-sequence stars we define the selection 
criteria separately for stars with T e g larger and smaller 
than 5250 K. For the stars with T e ff higher than 5250 K 
(displayed by the blue dots in Figure [2]), the selection 
criteria are as follows: 

• 5250<T eff <7000 K, 

• logy > (3.75-4.5)/(7000-5000)*(T eff -5000)+4.5, 

• (-Teff/500 + 14.1)< Mj <(-T eff /550 + 15.0), 
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where Mj is the estimated absolute magnitude in J band 
from Carlin et al. For stars with T e ff<5250K (presented 
by the red points in Figure [2]), the criteria are as follows: 

• 4000<T e ff<5250 K, 

• log g > 4.2, 

• (-T eff /700 + 11.0)< Mj <(-T eff /750 + 12.0). 

The criteria on log g and Mj are empirically set to ensure 
that the selected stars are in the main-sequence. The cri¬ 
teria for higher T e g remove the possible sub-giant branch 
stars so that the rest of the warm main-sequence stars 
may better trace the young populations. The criteria for 
lower T eS exclude the giant stars. We further select stars 
whose spectra have signal-to-noise ratio higher than 10 
in the g band, which is the most sensitive portion of the 
spectrum to the stellar parameters. We also select the 
stars within 

• 100< \z\ <500pc, 

• 7< Rgg <9kpc, 

adopting the position of the Sun to be (X, Y, Z) = (8, 
0, 0) kpc. The stars within 100 pc are not well ob¬ 
served in the LAMOST survey because most of these 
stars are too bright. Therefore, we only select the stars 
with |z| > 100pc. Finally, a total of 209,316 FGK 
type main-sequence stars are selected after applying all 
above criteria. It is worthy to point out that this sam¬ 
ple is significantl y extended compare to the previous 
works (DB 98 : S12:lNord strom et al.ll2004 iSiebert et al. I 
120081 : iBiidenbender. van de Yen fc Watkinsl 12014 etc.). 
Therefore, the statistics based on such a large sample 
may be substantially improved and the sample contain¬ 
ing a wide range of effective temperatures allows us to 
map the detailed kinematics for variant spectral types. 
Figure [3] shows the distribution of the sample in the 
Cartesian coordinate system with respect to the Galac¬ 
tic center (the top-left, top-right and bottom-left panels) 
and the T e g vs. \z\ plane (the bottom-right panel). The 
five vertical dashed lines in the bottom-right panel di¬ 
vide the effective temperature into six bins, in which the 
stellar kinematics are determined separately in the later 
sections. Figured] presents the distribution of the sample 
in l vs. b plane. From these figures, it is clear that the 
stars are not uniformly distributed in the sky. This leads 
to a distortion of the measured velocity ellipsoid, which 
is discussed in the Section 13.3.21 Figure [5] displays that 
most of the selected stars are with the signal to noise 
ratio in g band larger than 10. 

2.3. The radial velocity 

The LAMOST catalog provides the radial velocity of 
each star with uncertainty of about 5 km s -1 dGao et al. I 
120141) . In order to investigate the systematics of the ra¬ 
dial velocities, we cross identify the LAMOST data with 
the A POGEE data released in SDSS DR10 (jAhn et al. I 
120141 1 and find ~1000 common stars with good parameter 
estimates and velocity scatter smaller than 0.3 kms^ 1 
(which removes most multiple or variable stars) from 
APOGEE and signal-to-noise ratio higher than 20 in 
the LAMOST spectra. The radial velocity derived from 
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Fig. 2.— The distribution of the LAMOST DR1 stars in T e ff vs. 
log g plane. The black points are all stars in the LAMOST DR1 
with signal to noise ratio of SNR > 10 in g band . The blue and 
red points are the main-sequence stars selected for this work. 
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Fig. 3.— The distribution of the sample in Galactocentric 
Cartesian coordinates. The top-left (x-1^|), top-right (y-|z|), and 
bottom-left (x-y) panels shows the projected distributions, respec¬ 
tively. The Sun is located at (cc, y , z) = (8.0, 0, 0) kpc. The 
bottom-right panel shows the stellar distribution in Tea vs. |z| 
plane. The five vertical dashed lines divide T e g into six bins, in 
which the velocities and their dispersions are determined in this 
work. The contours represent the number of stars with 
different levels. 
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Fig. 4.— The distribution of the sample in l vs. b plane. 
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where N is the number of the stars. This can be specified 
as 

af oa = cr^ cos 2 Zcos 2 6 + Oy sin 2 1cos 2 b + cr^ sin 2 5 

+ a^r V sm2lcos 2 b + cr 2 w cos/sin26 (4) 

rr vw smlsm2b, 

where the dispersions <jjj , ay and aw specify the size and 
shape of the velocity ellipsoids, and the cross terms auv, 
ajjw and &vw determine the orientation of the velocity 
ellipsoids. 


Fig. 5.— The distribution of the signal to noise ratio in g band 

t lie'lliS/RfeT pipeline is slower by 5.7 kms^ 1 compared 
with APOGEE, as shown in Figure [G] The reason for 
this offset is unclear and is worthy of further investiga¬ 
tion in future works. In this work, we simply add an 
additional 5.7 kms -1 to the LAMOST derived radial ve¬ 
locity to match the other survey data. It is noted that 
the offset is weakly anti-correlated with T e ff, but 
no significant correlation with [Fe/H] and log g. 
We further investigate whether the weak anti¬ 
correlation with T e ff changes our result in sec¬ 
tion 15.41 

3. THE METHOD 

For most of the LAMOST dat a, the proper motions 
can be found from either UCA C4 (Zacharias et aTll2013 ) 
or PPMXL dRoser. Demleitner fe Schilbachl I2010T) cata¬ 
log. However, because the purpose of this work is to 
derive the velocity ellipsoids, it requires accurate mea¬ 
surement of the uncertainties of the proper motions as 
well as the values themselves. Unfortunately, after ten¬ 
tatively measuring the three dimensional veloci¬ 
ties and their ellipsoids with the combination of 
the proper motions and the radial velocity, we 
find that the errors of these proper motions are too large 
to derive the reliable velocity ellipsoids. Therefore, we 
turn to measure the three dimensional velocity and its 
ellipsoid from only the one dimensional line-of-sight ve¬ 
locities of the LAMOST sample spanning a large area of 
the sky. In this section, we first describe the velocity de¬ 
projection method and then discuss the validation and 
calibration of the approach. 

3.1. The de-projection method 

Consider a star with the line-of-sight unit vector r and 
heliocentric velocity of v. The observed line-of-sight ve¬ 
locity vios can be written as 

v ios = r v. (1) 

r is determined by the Galactic coordinates l and b and 
the three components of v in the Galactocentric Carte¬ 
sian coordinates are U, V, and W, respectively. Then 
Equation m can be expanded as 

vios = Ucoslcosb + Usin/cos6 + Wsinfc. (2) 

Given a group of stars at similar distance along the line 
of sight r, the velocity dispersion projected onto the line 
of sight should be 


3.2. The likelihood 

According to Equations [2] and U one is able to recon¬ 
struct the 3-D mean velocity and velocity ellipsoid at a 
given spatial position from the line-of-sight velocities of 
a group of stars located at the same position. To do this 
we construct a likelihood as the following 

(vios,k Vmodel^k) 1 
2 ( a model,k + £ l) _ J 

( 5 ) 

where vi os ,k is the measured line-of-sight velocity of fcth 
star, €ios,k is the measured error of the line-of-sight veloc¬ 
ity, v mod ei,k and a model,k are the predicted mean line-of- 
sight velocity and the projected velocity dispersion along 
the line-of-sight of the £;th star, respectively. The follow¬ 
ing 9 quantities, U, V, W, au, ay, aw, guv, guw, and 
ayw are free parameters in Equation ([5]). The best fit 
values for them can be found by maximizing the likeli¬ 
hood L. 

We use Markov chain Monte Carlo (MCMC) simula¬ 
tion to estimate the posterior probabilities of the veloc¬ 
ities and the velocity ellipsoids. The priors of the three 
mean velocities, U, V , and W, are evenly distributed 
in all real values, while the priors of the velocity ellip¬ 
soids are evenly distributed in the range of (0, + 00 ). In 
this work, we run the MCMC simulation with emcee, 
which implements t he affine-invariant ensemble sampler 
([Goodman fc Wearell2010f ). 

3.3. Validation of the method 

The velocity de-projection technique has been applied 
by DB98, who used the proper motions of Hipparcos 
to estimate the velocity and dispersions of stars with 
median dist ance of ^70 pc around the Sun. Later, 
iFuchs et al.l (120091 ) applied the similar method to the 
SDSS data, which extends t o a few hun dred p arsecs away 
from the Sun. However, McMillan & Binnevi (120091 ) then 
pointed out that the de-projection technique may pro¬ 
duce systematic bias in the cross terms, especially in 
ajjw, if a is a function of R and z. Because our method 
is very similar to that DB98 used and our data is located 
between 100 and 500pc in \z\, we may also suffer from 
similar biases. Furthermore, unlike the Hipparcos data, 
the highly uneven spatial distribution of the LAMOST 
data may lead to another systematic bias. Therefore, 
before applying this method to the LAMOST data, we 
need to carefully validate it and understand these issues. 
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3.3.1. Validation with GCS data 


N 











































5 




^ 10 

i 



logg, LAMOST (K) 


Fig. 6. Left panel: The offset of the radial velocity between LAMOST and APOGEE as the function of the LAMOST 
T e ff. The black dots are the individual stars and the red filled circles are the median values in each T e ff bin. The red line 
indicates the best linear fit of the median points. The black dotted-dash line marks the location of the average offset of 
5.7 kms -1 . Middle panel: The offset of the radial velocity as a function of the LAMOST [Fe/H]. The symbols are similar 
as the left panel. No clear correlation is found in this plot. Right panel: The offset of the radial velocity as a function 
of the LAMOST log g. The symbols are similar as the left panel. 


Our first validation uses data from the Geneva - 
Copenhagen Survey (GCS; iNordstrdm et al.1 l2004l h 
which is a volume-complete set of F- and G-type stars 
within ^40 pc, and provides three dimensional veloc¬ 
ity components for each individual star. This dataset 
does not extend to a large volume, hence it may not 
be severely affe cted by the systematic bias explained by 
[McMillan fc Binnevl ( 2009 ). Thus we expect that the de¬ 
projection method may give estimates close to the true 
values. 

A sub-sample of 3712 dwarf stars for the validation 
test is selected from GCS data with 5500 K<T e ff<6000 
K, \U\ < 200kms -1 , |Vj < 400 kms -1 , and \W\ < 
200kms -1 . Outliers beyond 4 a in the line-of-sight ve¬ 
locity distribution are removed. Only the line-of-sight 
velocity of these data is used to derive their mean ve¬ 
locities and velocity ellipsoid according to Equation 10 
via MCMC simulation. In order to reduce the effect of 
the uneven distribution on the sky, we separate the sky 
into small equal-area bins and randomly select roughly 
equal numbers of stars in each bin so that the samples 
are approximately evenly distributed on the whole sky. 
We repeat the arbitrary draw of the data in equal-area 
bins 40 times and estimate the velocity and its ellipsoid 
for each random drawing with MCMC. The final derived 
velocities and velocity dispersions are the averaged values 
over the 40 random draws. Uncertainties on the veloci¬ 
ties and their dispersions are composed of two parts: (i) 
the internal error from the MCMC; and (ii) the external 
error due to the selection of the sample. The total uncer¬ 
tainties of the derived values are the square root of the 
quadratic sum of the two parts. We compare the best fit 
values with the true values directly measured from the 
U, V, and W of the individual stars (see Table [T|)and 
find that the derived velocities and the ellipsoid are in 
good agreement with the true values. This confirms that 
the de-projection method can re-construct the 3-D kine¬ 
matics for the stars within a very local volume. 

3.3.2. Validation with the mock data 

LAMOST data are neither volume complete nor do 
they completely cover the whole 4-7r sky area. This sub¬ 
sequently leads to large spatial distortion in the veloc¬ 
ity de-projection method since many lines of sight are 
not observed. Moreover, the LAMOST data in this work 


spread from 100 to 500 pc in \z\, far beyond the volume of 
both GCS and Hipparcos data. The systematic bias due 
to the spatial variation of cr may be stronge r in this larger 
volume, as pointed out bv lMcMillan fc Binnevl (120091) . In 
order to investigate both sources of systematics, we intro¬ 
duce a second test, assigning mock line-of-sight velocities 
from pre-defined velocity ellipsoids to data with the same 
spatial position as actual LAMOST stars. The random 
draws of the line-of-sight velocities are based on a pre¬ 
sumed Gaussian distribution function in Galactocentric 
cylindrical coordinates with mean velocity components 
of [(uij), (vtj,), (u z )] = [ 0 , -26, 0 ] kms -1 and velocity dis¬ 
persions of 

< t(R,z ) = (cr Ro + 20z, <T 0 O +2C )z, cr Zo + 3CU) kms -1 

x exp[(i? 0 - R)/R a \, (6) 

where Rq = 8.0 kpc is the Galactocentric radius of 
the Sun, R a is set to 5kpc, cr 0 = (, <J Ro ,<j ( j, 0 ,a Zo ) = 
(45, 32,24) kms -1 following iMcMillan fc Binnevl (I2009T ). 
and 0 is in kpc. The cross terms t t r cr Rzi and <x| z are 
set to 0. For each observed star, its stellar parameters as 
well as the three dimensional positions are kept the same 
as the observed values, but its three dimensional velocity 
vector is randomly drawn from the Gaussian distribution 
function. The mock line-of-sight velocity for each star is 
then given by Equation 0, presuming that U = —v R , 
V = and W = v z in the solar neighborhood. 

The mock sample is first split out into three slices in 
height, i.e. 100 pc < \z\ < 300 pc, 200 pc < \z\ < 400 pc, 
and 300 pc < \z\ < 500 pc. The overlap between neigh¬ 
boring slices can help to smooth the results. Then each 
| z | slice is again split out into six bins in effective temper¬ 
ature with an interval of 500K ranging from 4000 K to 
7000 K. Here, we assume that the mean velocity compo¬ 
nents and velocity ellipsoids are the same for stars within 
the same T e s bin in each \z\ slice. The results are dis¬ 
played in Figure [7] with the solid lines representing the 
derived values and dashed lines for the true values. Be¬ 
cause the velocity dispersions are functions of \z\ and 
stars with different effective temperatures have different 
spatial distributions, the true values are not strictly flat, 
but show a slight bend with T e g. The reconstructed 
(U), (V ), and (W )are consistent with the true values. 
However, neither the velocity dispersions, nor the cross 
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TABLE 1 

Comparison of the velocity ellipsoids derived from different methods for the GCS data 



(U) 

(V) 

(W) 

Vu) 

Cv) 

(<rw) 

Cuv) 

{<?uw) 

(ervw) 

direct 

-9.86±0.64 

-23.49±0.46 

-7.574=0.33 

37.10±0.58 

25.69±0.67 

20.33i0.37 

11.04i0.97 

-6.37il.49 

4.05i2.39 

de-projection 

-9.73±3.01 

-22.84±0.74 

-6.67±1.22 

35.40il.02 

24.86i0.95 

18.97i0.84 

9.95i2.04 

-8.65i5.79 

3.17i6.88 


The first line are the results from the direct estimation, while the second row is the results derived from the de-projection method. 


terms are in agreem e nt wi th the true values. Although 
iMcMillan fe Binnevl (l2009f) showed similar deviation in 
the cross terms, their test did not demonstrate the bias 
in the velocity dispersions. The inhomogeneous distri¬ 
bution in sky position of the LAMOST data sample (see 
Figure 0, which can induce even stronger geometric dis¬ 
tortion, may be responsible for the deviation in the ve¬ 
locity dispersions. 

3.4. Calibration of the derived velocity ellipsoid 

In general, the distortion in the velocity ellipsoid due 
to the unevenly distributed sample can be modeled as 
skew and tilt on the true velocity ellipsoid. Therefore, 
the distortion can be corrected by a skew and tilt model, 
which can be mathematically presented as matrix trans¬ 
formations. 

The outline of the calibration is as follows. First, the 
mock stars are selected to be exactly located in the same 
spatial positions of the observed stars, combined with 
the simulated velocities based on a known velocity dis¬ 
tribution function, which is only qualitatively similar to 
the true one. Because the mock data are selected with 
the identical spatial distribution as the observed data, 
the true velocity ellipsoid of the mock data would be 
affected by the same skew and tilt as the observations. 
Thus, we can measure the extents of the skew and tilt 
from the quantitive differences between the de-projection 
derived velocity ellipsoids and their true values for the 
mock data. 

In the rest of this section, we denote that all quantities 
with subscripts - 0 , - m , -t, and - c represent for the observed 
values, the values measured from the mock data, the true 
values, and the corrected values, respectively. First, we 
correct the velocity vector with a simple offset: 

v c = v 0 —(Av), (7) 

where ( Av) = (v m — v t ), which is derived by averaging 
over six simulations. Table [2] lists the mean offsets at dif¬ 
ferent T e g bins and |z| slices. Most of the corrections are 
less than 1 kms -1 , implying that the distortion due to 
the unevenly distributed sample may not severely affect 
the derived mean velocities from the de-projection. 

Next we consider the velocity ellipsoids. For simplicity, 
we assume that the stars are in a three dimensional Gaus¬ 
sian distribution. It is worth noting that this assumption 
is not physically true, because the asymmetric drift of the 
stars in the Gal actic disk skews t he a zimuthal velocity to 
the slower side (lBinnev fe Tremainell2008l h An approx¬ 
imation of a Gaussian distribution may derive a slower 
peak of the azimuthal velocity distribution. However, the 
derived V and ay based on the Gaussian approximation 
with the GCS data shown in Table [T] do not show obvi¬ 
ous bias from the true values, implying that for the real 
observed data in the solar neighborhood, the Gaussian 


approximation may not contribute significant offsets in 
either V or ay. Therefore, we can still use a three di¬ 
mensional Gaussian profile to approximate the velocity 
distribution. A detailed investigation of the effect of the 
Gaussian approximation in the azimuthal velocity and 
its dispersion is discussed in Section l5Jl 
We denote £ t and £ m as the true covariance matrices 
and the one derived with the mock data from the de¬ 
projection, respectively. Both can be rotated to align 
with the given axis by the following transform: 

D t =Vi l V t Vt 

D m =V~ 1 '£ m V m , (8) 

where D t and D rr are diagonal matrices and Vt and V m 
are the rotation transform matrices. The skew due to the 
inhomogeneity rescales the axis ratios and hence can be 
quantified by the division of the two diagonal matrices. 

D t = AD m , (9) 

where A is the scaled matrix. The tilt due to the inho¬ 

mogeneity then can be quantified by the following trans¬ 
form: 

S t = VtAV^ m V m Vt\ (10) 

We define the correction matrices as A = V tAV” 1 and 
B = VmVf 1 . Then, the observed velocity covariance 
E 0 can be calibrated by 

S c = AS 0 B. (11) 

The calibration matrices can be derived from the simula¬ 
tions with mock data combining the observed 3-D posi¬ 
tions with arbitrarily drawn velocities from given velocity 
ellipsoids. In practice, the random mock data may intro¬ 
duce some statistical fluctuations in A and B. There¬ 

fore, we run six simulations to generate six sets of A 
and B. The final calibrated velocity ellipsoid is aver¬ 
aged over the calibration by the six sets of transform 
matrices. Figure [5] shows the calibrated results, which 
are now consistent with their true values after the cor¬ 
rection. It shows that not only the dispersions but also 
most of the cross terms have been well reconstructed by 
the calibration. The calibration corrects the systematics 
raised from both the unevenly spatial distribution of the 
observed sample and the spatial variation of a. 

4. THE RESULTS OF THE LAMOST DR1 DATA 

We apply the calibrated de-projection approach to the 
LAMOST data. Figure [9] and Table [8] show the results. 
The correlations of the three velocity and the six ve¬ 
locity ellipsoid components with T e g are demonstrated 
in three slices of \z\, i.e. 100pc < \z\ < 300pc (red), 
200pc < \z\ < 400pc (green), and 300pc < \z\ < 500pc 
(blue). The data at each \z\ slice are split out into 6 
T e ff bins between 4000<T e ff<7000K with an interval of 
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Fig. 7. — Comparison of the derived velocities and velocity dispersions with their true values for one of the six mock data sets. In all 
panels, the solid lines represent estimated values, while the dashed lines represent the true values of the mock data. The red, green, and 
blue indicate samples at 100 < \z\ < 300, 200 < \z\ < 400, and 300 < \z\ < 500pc, respectively. The left panels present the three velocity 
components ((£/), (V), and (W)) in Tiff bins in each \z\ slice. The second column panels show the three velocity dispersions and the right 
panels show the cross terms. 


TABLE 2 

The mean offsets of velocity to calibrate the geometric 

DISTORTION 


Teff 

(N) 

<A U) 

(AV> 

(AH') 



100< \z\ <300 

pc 

4250 

13229 

0.66±0.68 

0.224=0.89 

-0.094=0.41 

4750 

19593 

-0.104=0.25 

-0.234=0.29 

-0.234=0.34 

5250 

25805 

0.384=0.38 

-0.27±0.40 

-0.014=0.32 

5750 

33453 

-0.294=0.20 

0.154=0.42 

0.064=0.37 

6250 

23565 

0.364=0.31 

0.484=0.30 

-0.184=0.53 

6750 

6739 

0.244=0.79 

1.964=1.29 

-0.90±1.44 



200< | z | <400 

pc 

4250 

7600 

2.444=0.92 

-0.22±1.26 

-0.314:0.25 

4750 

16639 

-0.264=0.28 

0.234=0.73 

0.05±0.40 

5250 

27523 

0.014=0.32 

0.104=0.52 

0.01±0.23 

5750 

36303 

-0.03±0.27 

0.274=0.34 

-0.054=0.33 

6250 

18621 

-0.22±0.48 

0.634=0.63 

-0.064=0.46 

6750 

3436 

-0.774=0.92 

2.134=1.29 

-0.154=1.42 



CO 

o 

o 

A 

A 

o\ 

o 

o 

pc 

4250 

3797 

0.584=2.11 

0.044=1.17 

-0.30±0.48 

4750 

11443 

-0.164=0.43 

0.474=0.76 

-0.384=0.51 

5250 

22503 

-0.414=0.50 

0.044=0.59 

-0.01±0.26 

5750 

33595 

0.064=0.27 

0.204=0.40 

-0.404=0.30 

6250 

13986 

-0.244=0.91 

0.584=0.44 

-0.734=0.29 

6750 

1566 

-1.814=2.24 

1.364=1.83 

0.444:2.09 


The (A U), (AV) and (AW) are the average offsets of the ve¬ 
locities for various T e ff bins at different \z\ slice. 


500 K. Similar to the procedure for GCS data described 
in section 13.3.11 the uncertainty of the velocity and its 
ellipsoid is the square root of the quadratic sum of two 
parts: (i) the internal error from MCMC; and (ii) the 
external error from the selection of the sample data on 
the sky. To measure the error from the second part, 
we run a 40-time bootstrap to randomly resample the 
observed data on the sky. For the mean velocities, the 
uncertainty is mainly from the first part, while for the 
second moments, the uncertainties are dominated by the 
second effect. In this sense, the final uncertainty of the 
velocity ellipsoid cannot be effectively reduced until the 
LAMOST survey has covered a large, contiguous region 
of the northern sky, as it will by the end of the mission. 

4.1. The mean velocities 

The left panels of Figure [9] present the results of the 
three components of the average velocity (([/), (V), and 
{W)). 

The most obvious feature is that the three mean ve¬ 
locities shown in the left panels are all correlated with 
effective temperature for all \z\ bins. For the radial ve¬ 
locity (U) in the top-left panel, at X , e ff=4250K, (U) is 
around 10 — 12 kms” , which is at around 0 kms” 
with respect to the local standard of rest (LSR, the 
black dashed line in the panel) adopting the solar mo¬ 
tion of U® = 11.1 kms - (iSclionrich. Binnev fe Dehnen I 
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Fig. 8.— Same as Figure 0 but the estimated values have been calibrated via the method outlined in Section 13.41 


12010T> . In other T e g bins, (U) for the stars within 
100 < \z\ < 300 pc are approximately flat at around 
—9.07 to —8.74 kms" 1 . For the stars in 200 < \z\ < 400 
(300 < \z\ < 500), on the other hand, (U) increases from 
-8.98 ± 1.47 (-9.37 ± 2.04) kms -1 at T eS = 4750K to 
-6.49 ± 1.54 (-7.24 ±2.09) kms -1 at T eS = 6750K. Al¬ 
though the change of (U) from cool to warm is small, 
the increasing trend is quite clear, especially when T e g> 
6000 K. The direction of (U) is toward the Galactic cen¬ 
ter; thus the top-left panel of Figure |9] indicates that 
either the warmer stars tend to move inward the Galaxy 
or the cooler stars tend to move outward. 

For the azimuthal velocity ( V) in the middle-left panel, 
the stars show a clear break at T e g= 6000 K. The stars 
cooler than 6000K are located at around -30 kms -1 , 
while those warmer than 6000 K go abruptly up to 
around -20 kms -1 . Thi s is an other version of Pare- 
nago’s discontinuity (iParenago 1119501 ) and has been re¬ 
discovered by DB98 at B — V = 0.61, whic h correspo nds 
to Fp ff~ 60001< according to iSekiguchi fc Fukugital 
( 2000 ). The cooler stars are on average older than the 
warmer stars and hence most of them have had a longer 
time to experience scattering, which can increase the 
velocity dispersio n and therefore in crease th e asy mmet- 
ric drift (|.Ienkinslll992t iBinnev fc Tremaind 120081 1 . This 
leads to a slower azimuthal velocity for the older (cooler) 
stars. Another interesting feature is that the LAMOST 
data shows an unusual gradient along \z\ for warmer 
stars. For the stars cooler than 6000 K, the azimuthal 
velocity decreases with \z\, while it turns out to increase 


TABLE 3 

The solar motion with respect to the LSR. 



kms 1 

U 0 

9.58±2.39 

Fo 

10.52±1.96 (hn=2.6 kpc, L=5.2kpc) 
10.05±1.98 (±R=2.5kpc, L=5.0kpc) 

IF© 

13.09il.85 (/i fl =3.0kpc, L=6.0kpc) 
7.01T1.67 


with \z\ for the stars warmer than 6000 K. The former 
trend is natural since the stars at higher \z\ need more 
vertical energy and hence slightly lose their angular mo¬ 
menta and rotate slower, given that the stars are in equi¬ 
librium. However, it is difficult to understand why the 
warmer, or younger, stars show the trend turning around. 

As did the other two velocity components, the vertical 
velocity (W) also shows a correlation with effective tem¬ 
perature. Compared to the zero velocity with respect to 
the LSR (shown as the black dashed horizontal line at 
7.25 kms -1 (|Schdnrich. Binnev fc Dehnen I [2010 )). the 
stars with T e g< 5500 K roughly stay around the zero 
point. However, the warmer stars with T e ff> 5500 K 
show higher ( W) than the zero point, i.e., they are 
moving toward the Galactic north pole up to (W)~ 
—4 kms -1 , which is equivalent to ±3 kms -1 with re¬ 
spect to LSR. The stars located at larger \z\ seem to have 
a smaller (W) than those located in lower \z\, although 
the trend is not quite statistically significant. 

The stars cooler than 6000K over all \z\ bins are on av¬ 
erage old. Since they do not show a significant gradient 
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Fig. 9.— The mean velocities and velocity ellipsoids estimated from the LAMOST DR1 samples. The symbols are same as Figure 171 
The three black dash lines in the left panels indicate the zero points of the three components with respect to the LSR according to 
ISchonrich. Binnev fc Dehnen I 11201011 . 


in (V) it is safe to assume that they are in equilibrium. 
Then, they can be used to determine the solar motion 
with respect to LSR. The solar motion in radial and ver¬ 
tical directions can be derived from 


Ulsr = ~{U) 

W LSR =-(W). (12) 

Thus, we directly obtain Uq = 9.58 ± 2.39 kms -1 and 
Wq = 7.01 ± 1.67 kms -1 by averaging over the (U) and 
(W) listed in the 3rd and 5th columns of Table[5]for stars 
below 6000 K in all \z\ slices. 

In order to estimate Vq, the asymmetric drift is re¬ 
quired and a dynamical model ha s to be introduced. 
We st art from Equation (4.228) of lBinnev fe Tremaine! 
( 2008 1. which gives the approximation of the asymmet¬ 
ric drift from the circular speed, velocity dispersions, and 
the stellar density. 


v s -V e = v a 




(3 In (Wft) 
d Ini? 


R d(v R v z ) 

vk dz 


(13) 

where R is Galactocentric radius, z the height above the 
plane, v c the circular speed, and v the number density of 
stars, while a bar indicates the average value. Because 
the data are very close to the Galactic mid-plane, the 
cross term o f vrv z does n o t si g nifican tly vary with \z\ 
dBiidenbender. van de Yen fe Watkinsl 120141) and there¬ 
fore the last term of the equation is very small and can 


be negligible. Assume that the stellar density v and erp 
are exponential functions of R: 

v^exp(-^) (14) 

<7u ~ exp{-R ), (15) 

where Hr and L are scale length for the stellar den¬ 
sity and (Ju , respectively. Then Equation (4.228) can 
be rewritten as: 

”• 3 + + 1 - ^ l- < 16 > 

Therefore, the solar motion in V can be estimated from 
V Q = -((V)+v a ). (17) 


Here, we set h R = 2.6 kpc following jJuric_eWrI 


(2008 

) or hn = 2.5 kpc following Schonrich & Binnev 

(2012 

) and Sharma et al. 

(2014|1. According to 

van der Kruit & Freeman 

given an isothermal disk. 

201l{), L is equal to 2 Hr 
Although this relation has 


not been confirmed for the case of our Galaxy, it is 
a reasonable guess. As comparison, the different case 
of Hr = 3.0 kpc and L = 6.0 kpc is also consid¬ 
ered. The final results, including all three compo¬ 
nents, are listed in Table [3] Our estimates of the so- 
la_r_ motion in three di mens ion s are in agreement with 
ISchonrich. Binnev fe Dehnen I (12010D and iHuang et al.l 
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TABLE 4 

The vertical gradient of the velocity dispersions at eaci-i 

T eff BIN. 


T e a 

K 

dcru /d\z\ 
kms -1 kpc —1 

day /d\z\ 
kms -1 kpc —1 

daw/d\z\ 
km s -1 kpc -1 

4250 

-0.2±1.4 

54.54=3.4 

20.74=0.7 

4750 

10.4±1.7 

22.94=2.3 

18.94=1.0 

5250 

12.7=1=2.1 

14.34=0.4 

13.14=0.5 

5750 

14.14=2.9 

30.74=1.2 

24.64=1.1 

6250 

21.84=0.4 

29.44=7.0 

22.34=2.2 

6750 

23.04=3.4 

5.84=0.5 

28.44=8.3 


4.2. The velocity ellipsoids 

The second column panels of Figure [9] show the veloc¬ 
ity dispersions as functions of T e g and \z\, and the right 
panels show the derived cross terms. The velocity disper¬ 
sions show clear trends either along T e g or \z\. First, all 
three velocity dispersions show an abrupt drop at around 
T e g~ 6000K. This sudden drop has been seen in all \z\ 
bins. It is in agreement with the turn-around point in 
(V) shown in the middle-left panel in Figure [9] The de¬ 
creased velocity dispersions for stars with T e g> 6000 K 
again imply that the young stars are less affected by scat¬ 
tering and hence are kinematically cooler. 

Second, the velocity dispersions show vertical gradients 
in most of the T e g bins. Particularly, the vertical gradient 
of au is correlated with T e g , as shown in the second 
column of Table [I] Although the vertical gradients of 
cry and aw are mostly larger that those of au, they do 
not show clear trends along T e g (also see Table [4]). 

The ratios of the velocity dispersions shown in Fig¬ 
ure m indicate the shape of the velocity ellipsoids. 
Columns 12 and 13 in Table [8] also list the values of 
the two ratios. Both the ratios drop with |^| in most T e g- 
bins. The only exception is at T e s= 6750 K, in which the 
ratios seem not correlated with |z| (Figure HO])■ 

With derived cross terms, we can also estimate the 
vertex deviation and the tilt angle at each T e g bin and 
\z\ slice. The definition of the vertex deviation is 

1 2a 2 

l v = — arctan( 2 uv 2 ), (18) 

2 ajj — a v 

and the definition of the tilt angle is 

1 2rr 2 

a = -arctan( / ). (19) 

2 a v — a w 

The right panels of Figure [9] show the derived cross terms 
of the velocity ellipsoids. Their values are listed in Ta¬ 
ble [8] Figure [Tl] shows the estimated l v and a as func¬ 
tions of T e ff and |z|. 

We find that l v varies with \z\. It is around 0° 
for stars with T e ff< 6500 K at 100 < \z\ < 300 and 
200 < 1 2 1 < 400 pc and then it becomes positive at 
300 < \z\ < 500 pc. The vertex deviation for stars with 
T e ff> 6500 K, as an exception, changes with \z\ in the 
opposite sense, i.e., it decreases with \z\ from zero to 
—6.7°, although the uncertainty is very large. On the 
other hand, the tilt angle does not show cle ar trends with 
either T e g or \z\, which is different with iBinnev et al.1 
(|20l4 hereafter B14), who found the tilt angle changes 
for the hot dwa rfs com p ared to the giants or cool dwarfs 
using RAVE (iKordooatis et al.l 12013 1 data. The reason 


of the discrepancy is hard to be investigated, because of 
the large difference of the volumes between the LAMOST 
and RAVE data and the difference in methodologies. 

4.3. Comparison of the velocity ellipsoids with other 

works 

Although the volume of our sample at 100 < \z\ < 
300 pc does not exactly overlap with that of the sam¬ 
ple of DB98, which only extends to 70 pc, we will com¬ 
pare our sample to that of DB98 under the assumption 
that the velocity dispersions do not largely change from 
\z\ <100 to 300 pc. The color index B — V used by 
DB98 can be converte d into eff ective temperature ac¬ 
cording to lSekiguchi fc Fukugitai (l2000t) . Specifically, the 
data at around B — V = 0.412 (the color index bin 
4 in DB98) corresponds to T e g~ 6600 K, the data at 
around B — V = 0.525 (bin 6 in DB98) corresponds 
to T e ff~ 6200 K, the data at around B — V = 0.582 
(bin 7 in DB98) corresponds to T e g~ 5800 K, and the 
data at B — V = 0.719 (bin 9 in DB98) corresponds 
to T e ff~ 5200 K. Because DB98 did not directly show 
the dispersions along the U, V, and W directions, but 
showed the dispersions along the major axis of the veloc¬ 
ity ellipsoids and the ratio of the major axis to the middle 
and minor axis, we also rotate the velocity ellipsoid show¬ 
ing in Tableland find the longest, middle, and shortest 
axes of the derived ellipsoids to compare with DB98. Ta¬ 
ble [5] shows the comparison of the largest velocity disper¬ 
sion a\ and the two ratios of < 71/02 and 01/(73 between 
DB98 and this work. In most cases, the velocity disper¬ 
sions with roughly the same T e g are slightly larger in our 
data than those in DB98 by up to ~ 3 kins' 1 . This is 
reasonable since the height of the data in this work is 
slightly larger than those in DB98. 

It also shows that the ratio 01/02 derived in this 
work is on average consistent with DB98, while 01/03 
is slightly smaller than DB98. In other words, compared 
with DB98, the velocity ellipsoid in this work is slightly 
broader in the z direction, which is probably because 
the LAMOST data extend over a larger range in \z\. 
Notice that the errors in axis ratios are larger in this 
work. These errors may be overestimated because they 
are propagated from the large errors in the cross terms 
via the rotation transforms, rather than deriving directly 
from the observed data. 

Surprisingly, for the stars located at 100 < \z\ < 
300 pc, the vertex deviation is around zero for all T e g 
bins. This is inconsistent with the results in DB98. 
First, it is noted that Sections |3.3.1l and l3.3.2l have shown 
that the cross terms of the velocity ellipsoids are the 
ea siest af fected in the veloci ty de-projection method, 
as iMcMillan fc Binnevl (120091 1 have stated. Both the 
method in DB98 and in this work are based on the 
velocity de-projection. Therefore, both the results of 
the vertex deviation are possibly affected by systematic 
bias. Because DB98 de-projected the velocity compo¬ 
nents from the tangental velocities measured from proper 
motions, while we derive them from the line-of-sight ve¬ 
locities, the systematic bias in both works may go to¬ 
ward different directions. Second, even though we have 
calibrated the velocity ellipsoid, the corrected auw may 
still slightly bias from their true values (see the top-right 
panels of Figure [HJ. As a consequence, it is quite hard 
to compare the vertex deviation between this work and 
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TABLE 5 

Comparison of the velocity dispersions between DB98 and this work. 


T e S 

B-V 

DB98 

This work (100 < \z\ < 300) 



V1 

04/0-2 

0-1 /0-3 

lv 

o-i 

o-i/<r 2 

o-i/o- 3 

ly 

K 

mag 

kms -1 



degree 

km s~ 1 



degree 

5250 

0.719 

*17 2n +i - 4i 

0 / .zu_ 0 93 

1 44+ u,i2 

2 f)4 +u ' 0i 

i 3 .it?:S 

40.81i5.06 

1.51i0.29 

1.76i2.05 

-1.63i5.90 

5750 

0.582 

37.64 till 

1 ni+°- 07 
1,0 — 0.18 

1 7 o+ 0.48 
-*-• ‘ °- 0.04 

I0.2±|;g 

39.97i5.11 

1.50i0.33 

1 .86il.21 

-1.80i8.98 

6250 

0.525 

32 . 931 J;?! \ 

1 c-i+ 0.13 
i -Oi_ 0 .i2 

2 1 Q + 0 - 64 

z - iy - 0.19 

1 q+6.0 
J " y —6.1 

32.97i5.77 

1.62i0.40 

1.91i3.11 

0.45i6.62 

6750 

0.412 

26.261°;®° 

1 ««+0.12 

i-seio.is 

2 in+ 0.52 
z - iO - 0.15 

10.21|;° 

26.05i5.47 

1.32i0.44 

1.56i2.92 

-0.27i9.63 


TABLE 6 

Comparison of the velocity dispersions between B14 and this work. 


Height 

T e S 

B14 

This work 



Ol 

<Tl/(T2 

cri/0-3 

01 

a\l<J2 

04 / 0-3 

pc 

K 

km s 1 

kms 1 

100 < \z\ < 300 

>6000 

30.47i0.70 

1.52i0.11 

1.95i0.08 

29.51i4.89 

1.47i0.21 

1.73i0.25 

<6000 

38.42i2.30 

1.52i0.05 

1.71i0.05 

40.71i0.54 

1.51i0.01 

1.80i0.05 

200 < h| < 400 

>6000 

32.70il.54 

1.46i0.19 

1.80i0.08 

31.20i5.53 

1.51i0.25 

1.67i0.09 

<6000 

43.49i2.77 

1.57i0.04 

1.75i0.01 

42.14i0.80 

1.43i0.09 

1.70i0.03 

300 < h| < 500 

>6000 

36.57i0.70 

1.40i0.18 

1.69i0.03 

34.23i4.78 

1.42i0.06 

1.61i0.16 

<6000 

48.76i2.50 

1.63i0.04 

1.67i0.07 

42.83i0.76 

1.32i0.15 

1.62i0.05 



Fig. 10.— The axis ratios of the ellipsoids in various temperature 
bins at each \z\ slice. 

DB98. This issue cannot be easily solved until the 3 
dimensional velocities of stars are provided. 

B14 analysed the kinematics of ~400 000 RAVE stars, 
and decomposed the sample into hot and cold dwarfs, 
red-clump and non-clump giants. For each of these 
classes, B14 provided an analytic model for the velocity 
ellipsoid at each point in the (R, z) plane. With the help 
of this model, we compared the velocity dispersions be¬ 
tween B14 and this work. Table |B] shows the comparison 
of the largest velocity dispersion cri and the two ratios 
of ay/oy and oy/cr 3 between B14 and this work in the 
different layers. The oy , ay/ay, and ay/ay between the 
two works are in good agreement with each other within 
la for most cases. However, for stars with T e g> 6000 K 
and 200 < \z\ < 400 and those with T e g< 6000 and 



4500 5500 6500 



Fig. 11.— The vertex deviation l v and tilt angle a in various 
effective temperature bins at each \z\ slice. 

300 < \z\ < 500, the two works have differences by about 
2 -a. 

It is not easy to compare the dispersions with S12 since 
they separate the data into metallicity bins rather than 
temperature or color index. Moreover, most of the data 
in S12 are at heights larger than 500 pc, which does not 
overlap with our data. However, a qualitative compari¬ 
son is still worthwhile since SI2 directly used three ve¬ 
locity components to derive the velocity ellipsoids. First, 
we compare the velocity dispersions at 300 < \z\ < 
500 pc and T e g= 6250 K (see Table [5]) with the data at 
—0.5 < [Fe/H] < 0.2 and (z) = —0.69kpc in S12, which 
have an = 38.8ilg kms -1 , oy, = 27.7lg4 kms -1 , and 
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TABLE 7 

The test results of the non-Gaussianity of V 


°d 

(V) 

a v 

F® a 

kms 1 

80 

-16.6 ±2.5 

12.7 ±0.5 

14.0 ±2.4 

110 

-20.2 ±2.7 

19.9 ± 1.1 

15.2 ±2.6 

140 

-32.0 ±14.1 

33.3 ±5.6 

22.7 ± 13.1 

°d 

(V) n 

°V,G 

W.G b 

kms 1 

80 

-15.0 ±2.3 

11.9 ±0.4 

12.5 ±2.3 

110 

-16.4 ±2.2 

16.7 ±0.5 

11.6 ±2.2 

140 

-22.2 ± 10.9 

23.5 ±4.4 

14.2 ± 10.3 


a: The mean and standard deviation of V is calculated directly 
from their definitions. The V© is calculated from the direct 
calculated standard deviation. No Gaussianity is assumed 
in these derived values. 

b: The mean and sigma values are calculated from fitting a Gaus¬ 
sian profile. And Vq,g is calculated based on the Gaussian 
sigma crv,G- 

a z = 22.4lg'y kms -1 , and find that they agree with each 
other. Then, we compare l v = 9.60° ± 2.56° with their 
result, which is 8.3° ± 3.4°. These two quantities are 
again consistent with each other. We also compare the 
tilt angle a = —1.73° ± 0.76° with S12, which obtained 
—5.6° ±2.1°. Although the direction of the tilt angles are 
same, the values have about 2er difference. This is prob¬ 
abl y because the tilt_angle i ncreases very rap idly with 
\z\ (jBiidenbender. van de Yen fo Watkins! 12014H and the 
data in S12 are located in larger \z\ than those in this 
work. 

5. DISCUSSION 

5.1. The effect of the non-Gaussianity of V 

It is well known that the V distribution is not Gaussian 
but has a long tail to the slow side; hence, the Gaussian 
approximation in Equation d5j) may not be correct. It 
is worthwhile to investigate the systematic bias in the 
derived velocity ellipsoids due to the Gaussian approxi¬ 
mation. 

We adopt the analytical form of the distri bution func¬ 
tion defined by Equation (15) of iCuddeford &; Binnevl 
<1991 . namely, 

f{R,VR,v^) ~ exp( - 7 $*e yv */ v °) 

xexp- vl + 2 ( 20 ) 

where R is the Galactocentric radius, vr the radial ve¬ 
locity, Vtf, the azimuthal velocity, and y = 8 / Rd- The free 
parameters are Rd, ad, and v c . We predefine three sets 
of the free parameters and generate three sets of mock 
data via Monte Carlo simulations. Then we use Gaus¬ 
sian approximation to estimate the mean v^, and a In 
all three simulations, we set R = 8 kpc, v c = 220 kms - , 
Rd = 2.5 kpc. vr follows a Gaussian distribution with 
zero mean and sigma of ajj = adexp(—R/Rd)- We also 
add Vq = 12 kms -1 in the simulations. The only vary¬ 
ing parameter is ad, which is chosen as 80, 110 , and 
140 kms - , respectively, for each simulation. We then 
compare with the directly estimated mean and standard 
deviation of without assuming a Gaussian profile. 
The Vq is also calculated according to Equations m 
and m given L = 5 kpc. The comparison is listed in 


Table [Tj This shows that ( V ) is underestimated by a 
larger amount as ad increases. On the other hand, the 
velocity dispersion estimates from Gaussian fitting are 
also systematically underestimated by a few kms -1 . As 
a consequence, the derived Vq : g from the Gaussian fit¬ 
ting is also slightly smaller than the true value. 

Therefore, we can infer that the (V) and ay from 
Equation © may be systematically underestimated. 
The same is true of the estimated Vq. However, it is 
worth noting that the velocity distribution function (1201) 
is different than that of the real Galaxy. In reality, this 
underestimation may be very small and can be negligible. 
Indeed, the validation test with GCS data shows that the 
reconstructed V and ay are only smaller than the true 
value within 1 kms -1 , as shown in Table [1] Therefore, 
the Vq estimates may be only slightly smaller than the 
unknown true value. Going from low to high |z|, the 
disk becomes hotter and a,i increases. Then, the (V) is 
underestimated at higher \z\. However, this systematic 
underestimation likely occurs in all T e g bins and conse¬ 
quently cannot be the reason for the opposite trend along 
\z\ for the stars with T e g> 6000K. 

5.2. The age-T e ff relation 

We reveal that the kinematic features for the stars with 
T e ff> 6000 K are significantly different from those for the 
cooler stars. In this section, we first give a theoretical 
estimation of the age-T e g relation and then we discuss 
the possible reasons why the kinematic features are so 
different for the stars with T e s> 6000 K based on the 
age distribution in next section. 

In order to derive the ag e-T e ff relation, we first need to 
know the star formation history of the Milky Way. We 
adopt t he star formation history shown in Figure 1 of 
iSclionrich fc Binnev I (I2009T) . which contains a star for¬ 
mation peak at about HGyr age and a star formation 
rate decreasing with time (see the top panel of Figure fl2l). 
Second, we use synthetic isochrones dMarigo et al. Il2008ll 
to select the range of T e g for the main-sequence stars at 
given ag e. We also adop t the initial mass function de¬ 
rived by iChabrierl (120031) to assign weight for the main- 
sequence stars at different temperature. Finally, we set 
up the stellar distribution in age vs. T e g- plane, as shown 
in the bottom panel of Figure |T21 The blue thick line 
shows the mean age at different T e s. It shows that for the 
stars with 6500 <T e s< 7000 K, the mean age is between 
2 and 3 Gyr and there are almost no stars with age larger 
than 4 Gyr. For the stars with 6000 <T e ff< 6500 K, the 
average age varies between 3 and 6 Gyr. For the stars 
with T e ff< 6000 K, the average age is around 8 Gyr with 
a very broad range covering all ages. However, since the 
star formation rate has a peak at around 11 Gyr, the stars 
with age older than 8 Gyr dominate this region. The T e g 
of the abrupt change in age is perfectly consistent with 
that of the sudden change in the mean velocity and ve¬ 
locity dispersions shown in the left and middle panels of 
Figure HU This implies that the Parenago discontinuity 
occurs at around 6 Gyr. This val ue is a bout 3 Gyr earlier 
than the result by iQuillen fc Garnettl (|2001h . 

5.3. The asymmetric motion of the young stars 

Simila r asymmetric motion is also reported in other 
works. iCarlin et al. I (|2013i ) found both the vertical 
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Fig. 12.— The top panel sh ows the star formation history used by 
ISchonrich &: Binnev I 120011 ). The bottom panel shows the stellar 
distribution in T e g vs. age plane for the main-sequence stars based 
on the Cha brier IMF (Chabricr 2003) and the synthetic stellar 
isochrones IlMarigo et al. 1120081 1. The color codes the normalize 
stellar distribution and the blue thick line indicates the mean age 
as a function of T e fj. The dashed lines indicate the separation of 
the T e fj bins in this work. 


and radial velocities are not at around zero with the 
LAMOST data. Although the data analyzed by them 
are located in a larger volume than the data we use 
in this work, they are qualitatively consistent with the 
kine matic feat ures o f the young stars in our samples. 

I Williams et al.l (120131) also found similar trends in the so¬ 
lar vicinity with the red clump stars observed by RAVE 
(|Siebert et al. Il2011f) . Particularly, they found a faster 
Vrf, (see their Figure 10), a negative vr (see their Fig¬ 
ure 12), which is equivalent with a net positive U in 
the Cartesian coordinates, and slightly positive vz (see 
their Figure 13) at the position close to the Galactic 
nrid-plane. Although the quantitative comparison is 
qu ite diffi cult si nce the spectral types of the stars in 
iWilliams et al.l (120131 ) are completely different and the 
spatial sampling is also different, the orientations and 
the values of the offsets in three dimensional velocity are 
quite similar to our results. In general, red c lump stars 
are th ought to be quite young and IGirardi fe Salarisl 
(120011) argued that the peak of the age of nearby red 
clump stars is only around 1 Gyr. Then, it is not sur¬ 
prising that similar asymmetric motion is found in the 
red clump samples and the younger stars from our data 
set. 

There are possibly three channels to explain the un¬ 


usual asymmetric motions for the young stars. We dis¬ 
cuss these three scenarios here separately. 

First, if the young stars are still not completely relaxed, 
they may keep the peculiar motion of the molecular cloud 
in which they are born. Assume that the young stars are 
formed in a giant molecular cloud containing totally 10 3 - 
10 4 new formed stars within a scale of ~200pc, then 
the relaxation time for the group of stars is between 0.5- 
4 Gyr. As shown in Figure [12j for the youngest stars with 
T e fj> 6500 K, even the maximum age is only ~ 2 Gyr. 
Therefore, it is possible that these young stars are still 
not completely relaxed and still display some kinematic 
features of their birth place. 

Second, it is likely that the young stars are perturbed 
by the non-axisymmetric structure in the disk, e.g., the 
central rotating bar, spiral arms etc. Some studies have 
argued that the perturbati on must affect both young and 
old stars (e.g. lFamaev et al~ll2005l ). However, in general, 
orbits of stars in nearly circular orbits and very close to 
the Galactic mid-plane can be approximated by three 
harmonic oscillations in azimuthal, radial, and vertical 
orientations. Subsequently, the relatively fixed frequen¬ 
cies of the oscillations make the stars more easily respond 
to the resonance induced by the rotating bar or spiral 
arms. Since young stars are mostly in near-circular or¬ 
bits and concentrated in the mid-plane, they are more 
easily perturbed by non-axisymmetric structures than 
the old stars. However, the bar and spiral structures 
are mostly in the Galactic nrid-plane and thus may not 
affect the vertical motion of the yo ung stars (though see 
iDebat tista 20141: iFaure et all 120141 for some mechanisms 
by which the bar and spiral arms can vertically perturb 
the disk). Therefore, it is hard to explain why the young 
stars also show a net positive value in W. 

Finally, the disk may h ave been pert u rbed by a 
merging event. Indeed , IWidrow et al. I (120121) and 
lYannv fe Gardneil (120131) found that the stellar verti¬ 
cal density shows wave-like features, which coul d be the 
densi ty wave excited by a minor merger dGomez et al. I 

Uni). 

Certainly, none of the scenarios discussed above can 
be easily ruled out. Further observational and theoreti¬ 
cal works are required to better address this challenging 
issue. 

An interesting question raised here is why previous 
studies of local kinematics based on the Hipparcos data 
do not show the asymmetric motion. It is quite hard to 
answer it. The different volume of the Hipparcos may 
be one reason. More detailed investigations with the fu¬ 
ture Gaia data including accurate distance estim ates as 
well as the 3-D velocities dPerrvman et al.ll2001l) may be 
needed to complete the picture of local disk kinematics. 

5.4. The impact of systematic bias in the distance and 
radial velocity 

Although Figure |T| does not show any depen- 
dance of the distance estimates on T e ff, some ex¬ 
treme tests, e.g., the distance is only over- or 
under-estimated for the warmer stars but un¬ 
changed for the cooler ones, are very helpful to 
verify whether the kinematic results are robust, 
The left and middle panels of Figure [T2] show the mean 
velocities derived from the distances with artificial 20% 
over- and under-estimation only for the stars with 
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T e g> 6000 K, respectively. Essentially, the asymmetric 
motions of the young stars still exist in both (U) and 
(W) with slightly shifted values. And the Parenago’s 
discontinuity and turnover phenomenon in (V) are still 
substantial. We also do not find any systematic bias due 
to the over- or under-estimation of the distance in the 
velocity ellipsoids. 

Figure [6] presents a weak anti-correlation be¬ 
tween the offset of the radial velocities and T e ff. 
This drives us to further investigate whether the results 
are stable when the offset varies with T e ff. We then 
fit the offset of the radial velocity with a linear 
function of T e g and correct the radial velocities 
for the individual stars with the T e g- independent 
offset. We find that this leads to slightly shifts in (U), 
(V ), and (W), as shown in the right panels of Figure [13] 
Although the difference in (U) between the young 
and old stars is vanished, The differences in (V) 
and (W) are even stronger. It may also shift the de¬ 
rived local motion of the Sun. For U®, it may decrease by 
~ 1 kms - , and for W®, it may increase by ~ 1 kms -1 . 
These shifts are within 1 -a un certainties a c cordi ng to Ta¬ 
ble [3] We also notice that Gao et al. (|2015h found 
that the offset is 3.8 kms - compared the LAM¬ 
OST data with PASTEL catalog (|Soubiran et all 
12010 ). implying that the calibration of the re¬ 
dial velocity may also depends on different cal¬ 
ibrators. We did another test with the offset of 
3.8 kms -1 and found the mean velocities shift by 
~ 1 kms -1 for all stars and the systematic differ¬ 
ences between the young and old stars are un¬ 
changed. Therefore, in order to keep simplifica¬ 
tion, we adopt the correction of the radial veloc¬ 
ity with a constant 5.7 kms -1 for all stars in this 
work. According to the tests with various cor¬ 
rections, we conclude that the difference in (V) 
and (W) between the young and old stars are real, 
while the mild difference in ( U) shown in Figure [9] 
may be insignificant. 


6 . CONCLUSIONS 

We use the nearby FGK type main-sequence stars se¬ 
lected from LAMOST DR1 data to estimate the three 
dimensional velocities and velocity ellipsoids from line- 
of-sight velocities alone. It has been known that the ve¬ 
locity de-projection technique can introduce systematic 
bias due to the spatial variation of the velocity disper¬ 
sions. Moreover, we find that the uneven spatial sam¬ 
pling can also affect the velocity dispersion as well. In 
order to derive the corrected velocity ellipsoid from the 
de-projection method, we calibrate it using a set of trans¬ 
form matrices estimated from simulations. The calibra¬ 
tion works well in the simulations and then is applied to 
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the observed data. 

We associate the derived velocities and their ellipsoids 
with the effective temperature and \z\ and reveal that the 
asymmetric motions of the stars in the solar neighbor¬ 
hood reported by previous works are mainly seen among 
warm stars with T e s> 6000 K, which are very young with 
average age of less than 4Gyr. These young stars ro¬ 
tate faster in larger \z\ than in smaller \z\. Meanwhile, 
they move up toward the north Galactic pole by 
about 3 kms -1 and probably radially inward to 
the Galactic center by a few kms -1 . The nature of 
the asymmetric motion is still not clear. With the older 
(cool) stars, we give estimates of the solar motion with re¬ 
spect to the LSR. We obtain (U®, V®, Wq)=( 9.58±2.39, 
10.75il.96, 7.01A1.67) kms -1 . 

We also derive the velocity ellipsoids and find that the 
young stars have significantly smaller dispersions than 
the older stars. The vertical gradient of the velocity dis¬ 
persions is larger in ay and aw than in ajj. On the other 
hand, au shows clear correlation with T e g but the other 
two dispersions do not. We confirm that the Parenago 
discontinuity occurs at about 6 Gyr by comparing the 
velocity dispersions with a simple star formation model. 

The derived velocity ellipsoids in this work are essen¬ 
tially consistent with those in DB98, S12, and B14. The 
velocity de-projection method may still slightly suffer 
from the distortion in the derived cross terms even after 
calibration. Therefore, more accurate estimations about 
the orientation of the velocity ellipsoid should be done 
with three dimensional velocities. 
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Fig. 13. — The mean velocities with different systematically biased distance and radial velocity. The left panels are the results with 
over-estimated distance by a factor of 20%, the middle panels are the results with under-estimated distance by a factor of 20% only in 
the bins with >6000K and keep the distances in other bins unchanged. And the right panels are the results with the offset 
varying as an anti-correlation function of TeS presented in Figure [6] in the radial velocity distribution. 
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TABLE 8 

The derived velocities and the velocity ellipsoids 


Teff (N) (U) {V) (W) a u <r v <Jw <?uv 

K kms -1 


100< \z\ <300 


13229 

-10.414=1.90 

-26.09±1.48 

-7.264=0.96 

41.254=1.84 

27.164=2.11 

22.54±1.10 

-6.514=10.23 

19593 

-8.744=1.49 

-25.964=1.52 

-7.174=0.92 

40.774=1.30 

27.124=2.14 

23.434=0.91 

-5.964=9.15 

25805 

-8.764=1.41 

-26.534=1.68 

-6.354=1.23 

40.794=1.60 

27.104=1.55 

23.174=1.26 

-5.164=6.28 

33453 

-8.86±1.40 

-26.454=1.62 

-5.354=1.19 

39.864=1.36 

26.684=2.00 

21.804=1.04 

-5.294=9.62 

23565 

-8.924=1.50 

-21.444=0.93 

-3.644=1.16 

32.974=1.17 

20.21±1.99 

17.364=0.95 

2.314=6.00 

6739 

-9.074=1.37 

-23.174=1.85 

-4.574=1.75 

26.054=1.05 

19.684=1.97 

16.744=1.86 

-1.184=9.26 


200< | z | <400 


7600 

-12.224=2.19 

-28.474=2.00 

-7.88±1.30 

40.794=2.19 

31.534=2.86 

24.83±1.44 

-2.114=12.69 

16639 

-8.98±1.47 

-29.054=1.91 

-6.90±1.00 

42.344=2.07 

28.694=2.23 

25.62±1.09 

-9.784=11.28 

27523 

-8.204=1.87 

-28.684=2.05 

-6.774=1.98 

42.734=1.50 

28.404=2.43 

24.62±1.05 

3.654=4.46 

36303 

-8.084=1.63 

-27.974=1.57 

-5.734=1.09 

42.194=1.47 

30.134=1.89 

24.614=1.33 

8.794=9.63 

18621 

-7.394=1.41 

-19.894=1.29 

-4.014=0.95 

35.024=1.50 

20.924=2.14 

20.29±1.09 

8.124=3.20 

3436 

-6.494=1.54 

-20.474=1.77 

-5.284=1.40 

27.274=1.39 

20.42±1.77 

16.94±1.34 

-4.244=10.13 


300< | z | <500 


3797 

-10.584=3.25 

-35.844=6.44 

-7.494=2.78 

41.214=10.82 

38.064=5.93 

26.674=6.46 

9.624=17.82 

11443 

-9.374=2.04 

-32.454=2.79 

-6.334=1.06 

42.854=2.12 

31.704=2.44 

27.204=5.43 

5.964=9.46 

22503 

-9.034=2.08 

-30.244=2.10 

-6.984=2.03 

43.334=2.80 

29.964=2.79 

25.78±1.56 

8.684=8.89 

33595 

-8.92±1.62 

-29.114=1.35 

-5.854=1.39 

42.684=2.17 

32.814=2.67 

26.71±1.39 

14.134=7.54 

13986 

-8.40±1.78 

-18.994=1.88 

-4.644=1.05 

37.324=2.01 

26.094=2.84 

21.814=1.32 

11.134=5.17 

1566 

-7.244=2.09 

-17.394=1.81 

-6.624=1.91 

30.644=7.05 

20.844=7.12 

22.424=6.33 

-7.454=9.12 



-4.21±7.61 

-5.40±8.12 

-4.83±6.01 

-9.70±7.09 

-2.70±12.51 

1.794=3.47 


-5.00±15.45 

3.57±7.98 

3.47±8.18 

5.634=5.14 

-4.414=6.65 

1.054=7.94 


1.524=0.09 

1.504=0.09 

1.514=0.07 

1.494=0.08 

1.634=0.10 

1.324=0.11 


1.834=0.07 

1.744=0.05 

1.764=0.07 

1.834=0.06 

1.904=0.07 

1.564=0.12 


-0.084=7.41 

-2.174=8.12 

-1.634=5.90 

-1.804=8.98 

0.454=6.62 

-0.274=9.63 


-0.854=3.92 

-1.504=4.05 

-1.194=8.58 

-4.804=3.86 

-0.534=4.25 

0.464=6.22 


-11.214=9.35 

-9.714=3.75 

-4.404=5.92 

-7.394=4.78 

-5.664=6.71 

3.214=4.88 


-2.244=8.92 

2.564=11.48 

1.804=7.40 

-2.154=7.90 

-0.444=5.57 

-1.414=7.43 


1.294=0.11 

1.484=0.09 

1.504=0.09 

1.404=0.07 

1.674=0.11 

1.344=0.10 


1.644=0.08 

1.654=0.06 

1.744=0.06 

1.714=0.06 

1.734=0.07 

1.614=0.09 


-0.684=10.41 

-5.564=10.81 

0.764=5.18 

4.994=7.87 

4.734=6.51 

-3.164=12.35 


-6.754=4.42 

-4.714=3.95 

-0.914=4.56 

-2.664=5.12 

-2.254=4.76 

1.294=7.38 


-11.104=10.02 0.344=9.78 

-7.464=4.18 0.074=9.94 

-6.244=11.04 -4.104=9.81 

-4.914=2.45 -2.654=8.72 

-5.264=3.04 -3.204=7.56 

7.574=11.32 2.714=5.20 


1.084=0.31 1.554=0.36 

1.354=0.09 1.584=0.21 

1.454=0.11 1.684=0.09 

1.304=0.10 1.604=0.07 

1.434=0.12 1.714=0.08 

1.474=0.41 1.374=0.36 


27.204=14.35 -7.014=7.32 

1.804=13.30 -2.904=3.53 

4.474=8.44 -1.844=6.28 

14.104=8.03 -1.254=4.87 

9.604=8.41 -1.734=5.93 

-6.714=12.79 7.364=10.43 





















